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Abstract 

We consider general model with factorization between macroscopic 
flow and matter decay distribution. We show universality of be- 
havior of elliptic flow at small pt for identified final particles in the 
symmetric nucleus collision. At high pt we compare non-relativistic 
and relativistic models for boosted decay distribution. In the relativis- 
tic models with distribution having power-like tail we show existence 
of elliptic flow saturation. This means that the elliptic flow value V2 
tend to constant at high pt. We discuss the importance of determina- 
tion of elliptic flow orientation which can help us to compare different 
models. For example, we introduce potential expansion model which 
can reproduce azimuthal asymmetry but based on the assumption 
that there are only cold strong forces and no thermalization. This 
model show opposite sign of elliptic flow orientation in comparison 
with thermal model. We also show that at small pt sign of elliptic 
flow can change due to sign of second derivative of decay distribution. 
This fact tell us that we must more carefully link elliptic flow sign and 
properties of the model. 
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1 Introduction 



Azimuthal momentum space distributions is very interesting effect 
widely studying at RHIC. Tlie strong elliptic ffow generated in non- 
central collisions was founded in the set of experiments. Basis method 
is study Fourier expansion of azimuthal particle distributions ^ . For 
study elliptic flow we interesting V2 which is calculated from azimuthal 
average (cos(2(/>)), where angle (p defined relatively to reaction plain 
which is must be determined for each event and defined as plain based 
on the beam axis and line joining center of the nucleous. Next, we 
must perform average over all set of events. 

There are many phenomenological explanation of the elliptic flow 
phenomenon. The main goal of the theoretical studying is explain how 
transverse plain asymmetry of matter in the reaction zone is trans- 
formed into asymmetry of the momentum space distribution. Exis- 
tence of the azimuthal asymmetry is striking evidence of reinteractions 
in the reaction zone which enforce matter into asymmetric radial ex- 
pansion. The modern most desired explanation of low pt elliptic flow is 
existence of the Quark Gluon Plasma (QGP) and thermalization stage 
during nucleus collision. This allow us to use hydrodynamic with rela- 
tivistic equation of motion d^T^^ . But it is necessary to know pressure 
as function of energy density and chemical potential. Pressure can 
be founded from QGP equation of state which unknown, especially 
near phase transition point. Many numeric calculations performed 
in the longitudinal boost invariant ideal fluid model with equation of 
state of ideal gas: P ~ T^. In additional dissipative effects, non- 
equilibrium transport equations and other improvements can be used 
with hydrodynamic model. Unfortunaly, we can not reach full phys- 
ical understanding because too much uncontrolled phenomenological 
assumptions but so small experimental observable variables. 

Recent experimental data from RHIC for Au-|-Au collisions for 
identified particles [3] show saturation of f 2 in the region pt > 2GeV, 
Fig. n Hydrodynamic calculation can not successful reproduce V2 be- 
havior in the saturation region. The origin saturation is still unknown. 
At high Pi was proposed a mechanism of jet quenching - highpt pertru- 
bative gluons propagated through asymmetric reaction zone and lost 
its momentum interacting with dense matter 0. Also proposer the 
contribution to V2 from nucleus classical fields in the saturation regime 
In our calculation we show that the saturation phenomenon have 
mostly kinematic origin, but saturation values depend on details of 
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matter motion due to reinteractions. 

In this paper we try to separate contributions to V2 from macro- 
scopic flow of matter u{f3) and from effect of the local matter decay. 
And try to describe low pt and high pt of V2 behavior in unified lan- 
guage. In our model for final particles spectrum we assume factoriza- 
tion between macroscopic flow and matter decay distribution /. In 
general case, u{(3) can be hydrodynamic flow and / be a thermal dis- 
tribution of specified hardron at the freeze-out temperature. So we 
can say that we add local decay matter effects to the hydrodynamic 
model. On the other hand, u{(3) can be color string transverse motion 
due its repulsion and / be standard color string fusion spectrum. In 
same sense we generalize Blast Wave model (Hj and we not concern- 
ing on real sense of flow u{(3) and decay distribution /. The model 
does not cover jet quenching because pertrubative jets produced at 
the early stage of collision where macroscopic flow is still not formed. 



2 Elliptic flow at low and high pt 

Consider elliptic flow in cm. system at at midrapidity, y = 0. As- 
suming boost invariance as usual we may reduce problem to the two 
dimension transverse plain calculations. Let u{(3) = u{f3, ip) be a den- 
sity of expanding matter depend on transverse flow velocity (3. We 
always assume that exists some fimax < 1 and u{P) = if \(3\ > f3max- 
For flnal transverse particle distribution we can write 

F(f, = F(p, „ = = (1) 

where f{p) = ~ transverse fusion spectrum of small piece of 

expanding matter, f/sip) - modifled fusion spectrum of fixed particle 
type according to it transverse motion with /3 velocity (one can speak 
about "blue shift"). With aim to exclude flnal spectrum mixing be- 
tween particles of different species we consider only elliptic flow for 
identified particles. Factor Z{p) defined to satisfy azimuthal normal- 
ization 

27r 

y'F(p,0)#=i (2) 



Z{p) = J(^J u0)f0{p)d'(5^ d4> (3) 
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Figure 1: The minimum-bias (0-80% of the colhsion cross section) V2{pt) for 
irg, A + A and 0. 



We consider only collisions of equal nucleus such as Au+Au. Then 
the distribution u{P) have very important symmetry properties which 
follow from the basic properties of the nucleus collision geometry. 

u{P,-'tP)=u{(3,i; + 7T) = u{P,xlj) (4) 

This automatic gives 

/ u{(3,il!) sm{nip)d'ip = (5) 

JO 

/ u{p,ij)cos{{2n + l)ip)dip = (6) 
Jo 

So for u{f3) we have only even cos harmonics in the azimuthal Fourier 
decomposition. 

Elliptic flow can be calculated from Fourier decomposition of F{p, cf)) 
V2 = J F{p, cfy cos(2(/))# = (cos(20)) (7) 





4 



F{p,^) = ^{l + 2v2Cos{2cl)) + ...) (8) 

Now our task is to calculate V2 from ((TJ using 

The integrated V4 is about a factor of 10 smaller than U2 [2]- In 
cases when ui (directed flow) V4 and other higher harmonics enough 
small relatively to V2 (one can say that the azimuthal distribution have 
shape like cos 2<j)) we can simple estimate V2 with sufficient accuracy 

V2 = ^(^F{p,0)-F{p,^)^ (9) 

For find //^{p) in relativistic form we must boost lab p to the co- 
moving frame. Decompose p into (5 and (3 orthogonal directions 

p = pcos{(j)-i{;)-^ +psm{(j)-'ip)^ (10) 

Now we can easy perform boost in /? direction. In relativistic form 
for fisip) we have 



Up) = f 



/sin2(0-^) + A __ ^1 (11) 



We also need to relativistic transform of phase space d^p in 

d^p 



Performing calculation we obtain 



inv (12) 



T?i A.\ ^ f fa\f/^(^ Ppcos{(t>-ij)\ (3df3dip 

np.4>) = ^^jum,{p)[i- (13) 

where M - identified particle mass and distribution / must be con- 
sidered as function on p^ . 

Consider V2 behavior at small p. In this paper for simplicity we 
show calculation with //^(p) in non-relativistic form 

fp{p) = f{\p- Mpf) = f{p^ + M^(3^ - 2MPpcos{4> - V')) (14) 

Relativistic calculation based on (|llj) also can be performed. The 
resulting p'^ behavior will not changed. Only constant factor definition 
will changed. 
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Due to existence factor Z~^{p) in any (/)-constant multiplier 
in / u{p)f/3{p)cP (3 does not contribute to F{p, (p) and can be dropped. 
Distribution f/3{p) at small p can be expanded in Taylor series. Zero 
expansion order does not contribute to V2 because there is no depen- 
dence on (f). In zero order Z{p) does not vanish 

Z{p) = Zo = J u{P, xl,)f{M^[5^)[3d(3di, (15) 

In the first order we have / ti(V') cos((/> — = due to © and 

©. Second order have contribution to V2 

Fip,<l>) = Yg)j uiP,ij)\f"{M^P')i2PMpfcos\<l)-^)Pdpd^P 

(16) 

where we omit all terms which give zero contribution to V2- 
If we need only V2 then 

cos^(0 — ip) = cos^ (/>cos^ + sin^ (psin^ ^ ~ 2 "^^^ 2*?^^°^ ^V' (17) 
Performing all calculations we find 

V2 = C{{u},{f},M)p^ (18) 

where 

27rM^ J u{P, ip) cos 2iPf"{M'^p'^)P^d/3diP 

Ju{p,i;)f{M^P^)pdpdi> ^ ' 

We found that the V2 grow as at small p and vanish at p = 0. 
This fact is not depend on detail form of f{p) distribution. In the Fig. 
^we can see p^-grow of V2 for identified particles at small p. 

Consider sign of V2- In most cases / u{f3, ip) cos lipdtp > because 
we set zero azimuthal angle to flow maximum. For any known dis- 
tribution /"(p^) > for p > pq, where po - same specific sufficient 
small momentum. So we can conclude that in our model f 2 > as 
observed^. 

Now consider case when p is large. For reason of explicit calculation 
we assume that the f{p) at large p have power like tail and have 
asymptotic This is consistent with any known fits of particle 

distributions. At first we consider non-relativistic case. 

..J_/„<,V-(l.(f)^<f)co,,.-.,)"°.^. 

(20) 

In the low energy collision mean (3 can be small. So f"{M'^(3'^) and V2 may be negative. 



6 



We can use Mj3/p as small value and expand (^0]) in Taylor series. 



Z{p) 



cos 



V-) (21) 



where we omit first expansion order because it does not contribute to 
In leading order we can calculate Z{p) 



For cos^( 



Z{p) = 27rp-2" j u{f5)(f(5 (22) 

ip) in (|21|) we can apply rule (|17)) .Then we have 
'M(3a\'^ J u{(3) cos Itpd^p 



V2 



P 



2ju{(3)(Pl3 



(23) 



So we see 1 /p'^ decreasing of V2 in the non-relativistic model. 

In the relativistic model we must determine asymptotic properties 



of (ini) where / 

//3(p)=p-'" 



\P\ 



-2a 



0" 



1-0' 



■ COS^(0 — ■0) 



2/3 



1-/3^ 



cos((/) — V') 



(24) 



Inserting it into we see that the factor absorbed into 

factor Z{p) and does not influence on the V2- So V2 tend to some 
constant at high pt 



V2 



J u{(5)R (cos((/> - ^),(3) cos{2(i))(f (5d(t) 
J u0)R (cos((/> - V), /9) (P^dcj) 



(25) 



where i? (cos(</> — ^/^), /3) = f p{p) / p~'^°' ■ 

It is interesting to know real V2 behavior at high pt to distinguish 
relativistic and non-relativistic model. Unfortunaly, the current exper- 
imental data have large error bars at pt > AGev and high pt behavior 
in not clear. 



3 Elliptic flow orientation and directed 
flow 

Asymmetric nucleus collision is collision between to different nucleus 
at midrapidity or collision of equivalent nucleus at rapidity near zbyteam- 
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If collision is asymmetric then © for n = does not holds. And we 
can expect existence of the directed flow vi ^ 0. Non-zero vi is usu- 
ally linked with barion stopping, limited fragmentation and projectile 
nucleon bounce. 

Important property of directed flow is orientation relatively to re- 
action plain. Directed flow can lie only in the reaction plain due to 
symmetry. Of course, we assume that there is no symmetry breaking 
in each nucleus collision. If vi directions was determined, the elliptic 
flow f 2 must be calculated relatively this direction. If founded V2 posi- 
tive than we can conclude that the elliptic flow is in reaction plain and 
maximum matter flow in direction where almond shape reaction zone 
is thin. And vice versa, if V2 negative than elliptic flow in direction 
where reaction zone is thick. 

Consider potential expansion model. In this model we assume that 
there are only potential forces acting on produced matter after colli- 
sion. Each two near pieces of matter after its production have positive 
potential energy of reinteraction. Instead abstract "matter" we can 
think about strong color fields or about color strings. The origin of 
repulsive forces ab initio is not known. In pure phenomenological 
language we can address problem to non-pertrubative QCD vacuum 
properties. 

For simplicity we consider color strings placed in one dimensional 
line with equal spacing. Let point x = be a center of this string 
chain, p - string tension, u - repulsive potential energy of neighboring 
string per unit length. After time evolution of this system all potential 
energy transferred into kinetic energy of transverse string motion. Due 
to symmetry we can consider only strings at x > 0. Let index i 
enumerate strings from x = to positive direction. String having 
i = after expansion will not moving due to symmetry. Transverse 
rapidity 6*1 of i = 1 string 

u 

01 = arcch- - 1 (26) 
P 

Other strings rapidity is simple additive 

0i = i0i (27) 

So longest string chain have larger transverse velocity of leftmost and 
rightmost strings and in thick direction of reaction zone we have flow 
maximum. 
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Determination of u 2 sign relatively to vi orientation is very impor- 
tant for testing of validity of different models describing reinteractions 
during collision. For example, in the hydrodynamic-inspired models 
V2 positive due to maximum pressure gradient in thin direction. In 
the potential expansion model where role of matter play color strings 
we expect negative V2 due to maximum energy released in thick direc- 
tion. At different energies different mechanism can play a main role. 
First signal of V2 positivity at RHIC energies was found in 0. But 
this is not full story. In our model we can see in ()19() that V2 sign 
can changed not only due preferred expansion direction but due to 
negative sign of /". 

Near midrapidity at small y we can expand u{P, y) in Taylor series 

U0,y) = uo0) + ui0)y + ... (28) 

where uq{P) have zero directed flow and ui{P) - non-zero. 
In our model directed flow have form 

vi{p) = JzJf){J ^0)ff^iP)d'p) cos0# (29) 

Consider vi at low y and low p. Similarly (|18|) we get 

^^1 = Cyp (30) 

where C - some constant not depended on y and p. Note that the 
sign of vi have no sense due to impossibility of its measurement. Such 
y scaling we was observed in [2]. Dependence vi{p) is not measured 
yet. 

4 Conclusion 

In this paper we have studied azimuthal anisotropy at large and small 
Pt in generalized model with factorization between macroscopic flow 
and matter decay distribution. We showed universality of behavior 
of elliptic flow at small p for specified final particles - V2 grow as p"^ at 
small p and vanish at p = 0. This fact is not depend on detail form of 
f{p) distribution and matter flow n(/3). So we can say that it is model 
independent. 

At high Pt we compared non-relativistic and relativistic models 
of Lorentz boosted distribution /. In the relativistic models with 
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distribution / having power-like tail we show existence of elliptic flow 
saturation. In the non-relativistic model V2 decrease as More 
precise measurements of V2 at region pf > AGev are needed to make 
selection between models. 

Determination the sign of V2 relatively to directed flow is impor- 
tant criteria to for testing of validity of different models describing 
reinteractions during collision. To make conclusion about model se- 
lection we need experimental determination of V2 sign in wide range 
of collision energy. For example, we can compare hydro and potential 
expansion models. 

We also shown that at small pt sign of elliptic flow can change due 
to sign of second derivative of decay distribution /. This fact must be 
considered if we want to make model selection based on elliptic flow 
orientation. For example, for color string at sufficient small pt we have 
/" > 0. 

Acknowledgments 

We thank A.V. Dmitriev and N.V. Prikhod'ko for useful discussions. 
This work was supported by RFBR Grant RFBR-03-02-16157a and 
grant of Ministry for Education E02-3. 1-282 

References 

[1] S. Voloshin, Y. Zhang, Flow Study in Relativistic Nuclear Colli- 
sions by Fourier Expansion of Azimuthal Particle Distributions, 
arXiv:hep-ph/9407282 , 

[2] C. Adler et al. [STAR Collaboration], Azimuthal anisotropy at 
RHIC: the first and fourth harmonics, arXiv:hep-ph/0305265 . 

[3] J. Adams et al. [STAR Collaboration], Particle dependence of 
azimuthal anisotropy and nuclear modification of particle produc- 
tion at moderate px in Au+Au collisions at ^s^^^ = 200 GeV, 
,arXivm ucl-ex/0306007 

[4] Alex Krasnitz, Yasushi Nara, Raju Venugopalan. Elliptic Flow 
from Color Glass Condensate, arXiv:hep-ph/0209341 

[5] P. Huovinen, P.F. Kolb, U. Heinz, P.V. Ruuskanen and S.A. 
Voloshin. Radial and Elliptic Flow at RHIC: Further Predictions, 
arXiv:hep-ph/Q101136| 



10 



[6] Xin-Nian Wang. Jet Quenching and Azimuthal Anisotropy of 
Large pT Spectra in Non-central High-energy Heavy-ion Colli- 
sions, arXiv:nucl-th/0009019| 



11 



